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ABSTRACT 

We present some simple models to determine whether or not the accretion of cold dark 
matter by supermassive black holes is astrophysically important. Contrary to some claims in 
the literature, we show that supermassive black holes cannot significantly alter a power law 
density cusp via accretion, whether during mergers or in the steady state. 
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1 INTRODUCTION 

Over the past few decades, the idea of dark matter has become 
deeply rooted within the astrophysical community. In particular the 
ACDM (cold dark matter with a cosmological constant) theory has 
done very well in reproducing experimental results on large scales 
(see e.g. Bahcall et. al. 1999 or Sellwood & Kosowsky 2000). 

However, on small scales (galaxies and smaller) , the success 
of cold dark matter has been less striking ( Moore 200 1 ) . On galactic 
scales, current experiments which measure the velocity dispersion 
of a galaxy as a function of radius, have found enormous mass to 
light ratosjnaUj*atoae^rOTn^ & Borriello 

2000). While this is strong evidence for the existence of dark matter, 
the derived mass profiles do not compare wel l with those p redicted 
from simulations involving cold dark matter (Moore 2001 ). In par- 



ticular, cold dark matter simulations predict the presence of a central 
cusp in the density profil e, whereas the observed p rofiles seem to 
have a well defined core (Salucci & Borriello 2000|). This apparent 



misalignment of theory and experiment has lead to three main types 
of solution: 



1 . Dark matter is real but not cold ( Moore 200 1 ) 



2. Dark matter is not the solution and we should look to alter- 



native theories (Milgrom 198:) 



3. Cold dark matter is correct but current theories are missing some 
important physics. 

In other words: does the absence of a central dark matter cusp 
tell us something about the nature of dark matter, or does it tell 
us something about galaxy formation? Understanding this is es- 
sential for facili tating direct search e xperiments for dark matter 
such as DAMA ( Bernabei et. al. 2001 ). With any number of parti- 
cle candidates for dark matter and the difficulty of removing noise 
from direct detection experiments, particle physicists almost need 
to know what they are looking for before they can find it. As such, 



any complimentary method which can rule out some candidates is 
extremely important. Unfortunately, before any robust conclusions 
about the nature of dark matter can be made from the rotation curve 
measurements of nearby galaxies, we need to be sure that what we 
are measuring is the dark matter profile as it would have been before 
the formation of a galaxy. It is this profile which contains informa- 
tion about the detailed nature of dark matter. If the profile has been 
changed due t o astrophysical effects such as matter outflow (G nedin 



& Zhao 2001) or bars ( Athanassoula. 2001 ) then some or maybe 
all of the information about the nature of dark matter could have 
been lost. It is important to work out which astrophysical effects we 
should consider and which are not important. 

In this paper, we explore the fate of a CDM cusp in the presence 
of a central black hole. Recently, Zhao, Haehnelt & Rees 2001 
(hereafter ZHR) have suggested that tidal stirring from infalling 
satellites could lead to significant dark matter accretion, accounting 
for some 20-40% of the mass of a central supermassive black hole 
if loss cone refilling can be efficient enough. As such, we consider 
whether or not supermassive black holes can alter the dark matter 
distribution in a galaxy in any significant way. 

In section 2, we recall that the steady state accretion rate of cold 
dark matter onto a supermassive black hole is negligible. In section 
3, we present a toy model for calculating the mass of cold dark 
matter accreted during merger events. In section 4, we discuss the 
validity of our model and the subsequent effect on the underlying 
dark matter distribution. Finally, in section 5 we conclude. 



2 THE STEADY STATE ACCRETION RATE 

The steady state accretion rate of matter onto a black hole has been 
treated by several authors (e.g. Lightman & Shapiro 1977). Any par- 
ticles which lie on orbits which b ring them within the Schwarzschild 
radius of the hole (see equation |ll|) will be swallowed. Since both 
low energy, high angular momentum and high energy, low angular 
momentum particles can be swallowed the problem is necessarily 
two-dimensional. The fraction of particles in phase space which 
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will be swal lowed is described by the 'loss cone' (L ightman & 
Shapiro 1977) and any particles lying within the loss cone will be 
swept away in a dynamical time ~ the time taken for a particle to 
complete one orbit. 

As can be seen in figure [j], the loss cone can be thought of as 
the locus of points swept out from trying to point a particle at the 
black hole from some distance r away. 

Now, b y considering particle orbits in a Schwarzschild metric 
( Wald 1984 ), any particle with an angular momentum per unit mass 
less than L ma x = VTIGMbh /c will fall within the Schwarzschild 
radius and be swallowed (where Mbh is the mass of the black hole, 
c is the speed of light in vacuo and G is the gravitational constant). 

The solid angle subtended by the loss cone, £2 can be written 

as: 



sin 9'd8'd(j) 
n(v t /v) 2 = n(L max /rv) 2 



(1) 



So, given some phase space density of dark matter particles, 
f(r, v), the total mass inside the loss cone Mi c is given by: 



oo (T 2 oo 

Mic = I Q4'irr z 4irv 2 f(r,v)drdv 



I6n 3 Lmax 2 f(r, v)drdv 



(2) 



Notice that the lower limit for the integral over r is the 
Schwarzschild radius. This is because the black hole itself must 
provide a cut-off for the dark m atter distribution. The upper lim it is 
given by V2oo, the virial radius (Navarro, Frenk & White 199€ ). 

Assuming a Maxwell Boltzmann distribution for the velocity 
density of the dark mattej^articks^ini^ai^ 
Frenk & White 1996) for the spatial density gives us: 




Figure 1. Schematic diagram of a black hole and its loss cone. The solid 
circle at the origin represents a dark matter particle a distance r away from 
a black hole (depicted by the large circle). The dark matter particle has 
velocity, v and transverse velocity Vt = L max /r. 



the loss cone. We look at this idea in more detail in the following 
section. 



pNFw{r)g M B{v) 



= *! ( ? V e (^) (3) 

Where po and ro are the scale factors for a particular galaxy, 
and < v 2 > 1 / 2 is the velocity dispersion of the dark matter halo. The 
velocity dispersion is taken to be a constant, which is a reasonable 
approximation given the observed rotation curves in many galaxies. 

For a 10 6 Mq black hole at the centre of a Milky Way type 
galaxy (we take scaling values as in table |l|), the dark mass within 
the loss cone is ~0.6 Mq, while the total mass of dark matter 
is some 1.9 x 1O 12 M0. Thus a tiny fraction of the dark matter 
is accreted in a crossing time but then not much else can happen 
until the loss cone is refilled. For stars refilling can be achieved by 
two-body interactions, continuing satellite mergers and dynamical 
friction, allowing a continua l diffusion of matter into the loss cone 
(Lightman & Shapiro 1977). However, for cold dark matter, the 



interaction cross section is presumed to be so small that the time 
taken to re-fill the loss cone by diffusion is very much longer than 
the age of the universe. 

The problem can be avoided by postulating that dark matter can 
self-interact thus raising the interaction cross section a nd allowing 
significant steady-state accretion rates (Ostriker 2000). ZHR have 
also suggested that tidal stirring from infalling satellites could refill 



3 ACCRETION DURING MERGER EVENTS 

As shown above, the steady state accretion rate of dark matter is 
negligible. However, during a galaxy merger event, while a central 
supermassive black hole is stationary with respect to the dark matter 
in its host galaxy, it is not stationary with respect to the dark matter 
in the merger galaxy. In the rest frame of one of the galaxies, the 
central black hole, belonging to the other galaxy will spiral in, 
sweeping up dark matter on its way. As the two galaxies merge, 
t heir CDM h alos will be tidally stripped out to the Roche limit 
(Roche 18591 where the density of a halo approximately matches 



that of its surroundings. We can expect, then, that the CDM cusps 
from the two galaxies should make it right into the centre of the 
forming merger remnant. 

This idea is borne out by various numerical studies (see e.g. 
Van Albada & Van Gorkom 1977 & more recently Milosavljevic & 
Merritt 2002) which show that the initial conditions are statistically 
preserved in major merger events. As such, it is not clear that a dark 
matter cusp will survive during galaxy mergers. Two black holes 
could, in principle, mutually accrete dark matter from each other's 
host cusp. The magnitude of the effect will depend on the loss cone 
refilling rate during such a merger process and it is this that we will 
discuss in 3.1 and 3.2. 

The presence of a significant bulge component in galaxies, 
apparently correlated with the central black hole mass, can set a limit 
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on the efficiency of this process. If the cusp depletion is found to be 
very efficient then it must occur in the early stages of bulge growth, 
or in low surface brightness (LSB) galaxies where the dark matter 
dominates over the baryonic matter even in the central regions. 

Faber et. al 1997 examined the central density profile of 61 
elliptical galaxies and spiral bulges. They found that, while all of 
the galaxies in their sample were well fit by a power law profile, the 
brightest galaxies contained a central core in the light profile at a 
break radius of ~ 10 — 100 parsecs. They suggest, as have others 
(e.g. Quinlan 1996 and Milosavljevic & Merritt 2001, hereafter 
MM), that the central cusp depletion could be due to scattering 
during major merger events as two massive black holes form a 
hard binary. We consider a complimentary mechanism, operating 
earlier in the merger process, based on the accretion of scattered 
material during galaxy mergers. We show that this mechanism is 
not important for growing black holes or rearranging dark matter 
cusps. 

Although in our analysis we will talk mainly about the effect 
on the central dark matter, our analysis could be applied equally 
to the central baryonic matter cusps which are observed in bright 
galaxies. 

We treat the infalling black hole as a scattering body which 
spirals into the centre of a galaxy via dynamical friction. Some of 
the dark matter which it scatters will be scattered back onto itself - 
self scattering, and some will be scattered into the loss cone of the 
other hole - mutual scattering. The sum of these two effects will 
give us an estimate of the total accretion rate of dark matter. We 
now look at these effects quantitatively. 



3.1 Self Scattering 

Consider two in-spiralling black holes, A and B, from the refer- 
ence frame stationary relative to hole B. In the self scattering case 
we consider the in-spiralling black hole, A, as an effective cross 
sectional area sweeping through phase space. 

We need only consider the growth of hole A though, by sym- 
metry, hole B will also grow via this mechanism. We can then 
parameterise hole B solely by the nature of its host galaxy. This is 
because of the empirically observed correlation between bla ck hole 
mass and bulge velocity dispersion (Ferrarese et. al. 2001). Small 
black holes reside in shallow potential wells, larger black holes 
reside in larger potential wells. 

We assume that hole A is brought in via dynamical friction 
from the dark matter and baryonic matter in the centre of the sta- 
tionary galaxy (hole B's host galaxy). Following the Chandrasekhar 
prescription (c.f. Binney & Tremaine 1987), we have that the inward 



velocity of the hole, -j- is given by: 



dr 

~dt 



-KM A 



rp(r) 



(4) 



where, p(r) is the density distribution of the central matter 
(baryons and dark matter), Ma is the mass of hole A, v c (r) is the 
local circular velocity and K is approximately a constant. For a 
Ma xwell Boltzmann distribut ion of particle velocities, K is given 
by (Binney & Tremaine 19871: 



K 



A = 



0.428G 2 27rln(l + A 2 ) 



< tr > 



(5) 
(6) 



GM A 

where, G is the gravitational constant, b ma x is the maximum 



impact parameter for the in-spiralling hole and < v 2 > 2 is the mean 
velocity dispersion of the surrounding baryonic and dark matter. 

From equation Q we can derive the dynamical friction infall 
time tdyn for the hole starting at radius n. This gives us: 



"dyn 



1 v c (r) : 



dr 



(7) 



KM A rp(r) 

Equation ^ is the Chandrasekhar equation for dynamical fric- 
tion and is approximately valid provided that the infalling mass doe s 
not exceed the mass interior to its orbit ( |Binney & Tremaine 1987 ). 
Notice that the minus sign in equation |4j has been substituted for 
swapped integration limits. 

We can also use equation ^ to derive the amount of mass swept 
up by the black hole as it moves from a radius a to a radius b. The 
distance travelled by the hole in a time dt is given by: 



ds = dtv c (r) — 



v c (r) 4 dr 
Krp(r) 



(8) 



If the black hole cross section is given by a a (r) then the mass 
of dark matter swept up by the hole by moving through a distance 
ds is then: 



dM A = 



v c (r) 4 (TA{r)p(r)dr 
Krp(r) 



(9) 



and so the total mass swept up as the hole moves from a radius 
a to b is: 



AM A = / dM A = 



1 v c (r) 4 a A (r) 



KM A 



3.1.1 The black hole cross section 



dr 



(10) 



If the black hole cross section were just nrl, where r s is the 
Schwarzschild radius, then the volume swept up in phase space 
would be truly tiny. The Schwarzschild radius is given by: 



2GM A 



(11) 



where c is the speed of light in vacuo. So, for a 10 6 M Q black 
hole, r B ~ 10 -7 parsecs. 

However, a massive body such as a black hole interacts with the 
medium through which it moves, preferentially scattering particles 
onto orbits which will bring them within the Schwarzschild radius, 
allowing them to be swallowed. Taking this gravitational focusing 
into account gives us a total cross section of (B inney & Tremaine 
1987): 



2 , 2-KGM A r s 



« c (r) 2 

„2 



1 + 



v c (r) 5 



c?v c (r) 2 



(12) 



Since v c (r) must always be less than the speed of light, the 
cross section due to gravitational focusing will always be much 
more important than the physical cross section of the black hole. 

Putting this all together we get, for the mass swallowed as the 
hole moves from radius a to 6: 
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Figure 2. Radial dependence of the Coulomb factor, ln(l + A 2 ), in power 
law density profiles with exponent a. The solid line is for a = 0.5, the 
dotted line is for a = 1, the dashed line is for a = 1.5 and the dot-dashed 
line is for a = 2. 



AM A 
A 



AMA V -^ldr 



0.428 ln(l + A 2 )c 2 



(13) 
(14) 



Using equation 13 we can now calculate the mass swept up 
as the black hole spirals in. We cannot, however, naively integrate 
from some start radius to zero to obtain the mass swallowed. This is 
because, if the black hole manages to complete a whole orbit before 
moving inwards a distance greater than ~ a/ovl then we have the 
problem of shell crossing; that is, we would be counting the mass 
more than once. In order to avoid this problem, we can make the 
distance a to b some small interval and then numerically sum each 
mass shell Am, comparing each Am with the actual mass of cold 
dark matter present, and taking the smaller value. In this way, we 
avoid divergent mass consumption at small radii where the same 
part of configuration space is swept out many times but the mass 
within that space can be swept up only once. Quantitatively, this 
means that we require: 



AM A < 47T 



2 p(r)dr Va, b 



(15) 



3.1.2 The Density Profile 

The density profile for the dark and luminous matter in most galaxies 
can be well appr oximated, within a few kiloparsecs of the centre, 
by a power law (Blais-Ouellette & Carignan 2002). As such, we 



adopt the following form for p(r) 



P(r) = P o[- 



(16) 



For the dark matter, a is observed to be typically a < 1, 
whereas CDM simulations predic t dark matter slopes of a > 1 
(|Blais-Ouellette & Carignan 2002|). 



1987): 



The^rca^chnUarvdccit^i;^ & Tremaine 



Parameter 


Value 


Description 


Free Parameters 






M bh (M e ) 


10 6 


Black hole mass 


r,:(kpc) 


1 


Black hole start radius 


M 2OO (M ) 


1.9 x 10 12 


Halo mass within virial radius 


h 


1 


Hubble parameter 


Derived Parameters 






c„ 


13.34 


NFW halo concentration 


^2oo(kpc) 


201.4 


Virial radius 


p o (M kpc- 3 ) 
r' (kpc) 


2.49 x 10 7 


NFW scale density 


15.14 


NFW scale length 



Table 1. Initial Conditions. See text for details. 



v c (r) 



r' 2 p(r')dr' 



4nG 
r 



2 — a 



(17) 



< v > 



3 — a 

For a power law profile, it can be shown (see Appendix I) that: 

2 



(18) 



and if the maximum impact parameter, b ma x , is taken to be the 
starting radius of the hole, Vi, this gives us: 



4np ron 2-c 

A = T 

M A (3-a) 



(19) 



Figure | plots ln(l + A 2 ) against r for initial conditions set 
up as in table |l| and for varying a. For r > O.Olkpc, or a > 0.5, 
ln(l + A 2 ) ~ constant. However, in general it is not constant over 
the full range of r and hence neither is A (see equation 14). We 
assume that A is constant over the small change in r, Ar — a — b, 
but does vary as a function of a. 

Thus, putting this all together gives us, for the mass swept up 
by hole A: 



... 4TvGA(a)M A porS' . 2 - a 

AMa = ~~R \T^, ^ ( a 

(3-a)(2-a) 

AM A = 4TrGA(a)MAp r 2 ln(^) 



a < 3, a / 2 

a = 2 
(20) 



3.1.3 Initial Conditions 

From equation 20, we can see that the important factors for deter- 
mining the amount of mass accreted are Ma, the mass of the hole, 
and po, ro and a, the density parameters for the central cusp. 

The initial infall radius of the hole can be taken from numerical 
simulations performed by MM. We use lkpc as our starting position 
for the hole. 

The density parameters are scaled such that the galaxy has the 
same mass as an NFW-type profile within lkpc. The NFW profile 
is given by (Navarro, Frenk & White 1996): 



Pnfw 



lim pnfw 

r<<r' 



A) 



(^)(l + 7^) 2 



/ / 

r 



(21) 
(22) 



Where p' and r' are given by: 
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Pa = 

So = 

r = 

»~200 = 



Sopcrit 

200 c 3 n 
~ln(l+c n )- 

5*200 /c„ 

M200 
|7T200pcr 



(23) 
(24) 

(25) 
(26) 



Where noo is the virial radius, M200 is the mass within the 
virial radius and is taken to be the mass of the galaxy, c n is the 
NFW concentration parameter, and p cr it = 277.3/i 2 Mq kpc~ 3 is 
the critical density. 

The only free parameter in the NFW profile is Af2oo since the 
concentration param eter is found, numerically, to be a function of 
the mass of a galaxy (Navarro, Frenk & White 1996). 



The power law scale density and radius (po, ro) we use are then 
scaled to contain the same mass as the NFW profile within lkpc. 
Thus: 



/■o 



'•() 



Po 



ro 
lkpc 



(3 -a) 



(27) 
(28) 



Recall that this density profile describes the galaxy which is 
in the stationary frame. The infalling galaxy is parameterised only 
by the mass of its black hole. Thus, by setting up the profiles for 
a Milky Way type galaxy, we can model the full range of mergers, 
from equal mass to minor mergers, by scaling the infalling black 
hole mass. 

A summary of all of the initial conditions used is given in 
table |l|. The parameters are broken up into those which are put into 
the model and those which are subsequently derived. The scaling 
parameters for the NFW profile, for example, are set up consistently 
for a halo which contains a mass of 1 .9 x 10 12 Mq within the virial 
radiu s. A halo of this mass has a con centration parameter of about 
13.3 ( |Navarro, Frenk & White 1996| ). 



3.1.4 Results 

The numerical integration was performed using IDL, with the initial 
conditions set up as in table [l and with a numerical resolution of 
Ar = a — b = 0. 1 parsecs. The results are displayed in figure Q 

As can be seen, most of the accretion occurs in the outer regions 
and increases with a. However, even for a = 2.9, the effect is small, 
with less than IOOM0 accreted in total. This is a tiny fraction of the 
black hole's mass, which is some 10 6 M Q . For an NFW profile, with 
q = 1 (dotted line, figure ^), the total mass accreted is less than 
in the steady state case. Finally, from equation 20, we can see that 
the total mass accreted scales linearly with Ma, so the fractional 
increase in black hole mass is unaffected by any increase in Ma- 



3.2 Mutual Scattering 

Now we consider the case where the black holes scatter dark mat- 
ter into each others' loss cone. Consider, as in the self scattering 
scenario, a black hole, A, spiralling into a hole B, from a reference 
frame stationary with respect to the hole B. As before, each black 
hole has its own bound CDM cusp. The dark matter hole A scatters, 
for example, is that which resides in hole B's CDM cusp. 



T=1.46055e + 09 



10.0 
r(pc) 



Figure 3. Mass accreted by a 10 6 Mq black hole spiralling in to a Milky Way 
type galaxy from 1 kiloparsec away from the centre due to self scattering. 
The infall time, T, is marked in years. The solid curve is for a = 2.9, the 
dotted curve is for a = 1 and the dashed curve is for a = 0.5. 



Hole A will scatter some fraction, tja, of particles within its 
sphere of influence into the loss cone of hole B, and these particles 
will then be quickly swallowed in a dynamical time. Similarly, hole 
B will scatter particles into hole A's loss cone and hence both holes 
will grow via accretion as they merge. 

The problem, then, is similar to the self scattering case except 
that now, the effective black hole cross section is due to the sphere 
of influence of hole A, and hole A causes the growth of hole B. 

Since each hole induces the growth of the other, we should 
consider both Ma and Mb as functions of r and consider the 
growth of both holes. However, we will make the assumption that 

A , A 5 A ~ a * Jb ~ small so that Ma and Mb are ~ constant. 

ma Mb 

The sphere of influence of hole A is given by: 



Tsph,a 



GM A 

<v 2 > 



(29) 



1/2 

Where Ma is the mass of hole A and < v > 'is the local 
spherical velocity dispersion in hole B's host galaxy. For a 10 6 M Q 
black hole in a Milky Way type galaxy, r sp h,a ~ 0-3 parsecs. Thus 
the effective black hole cross sectional area, T^r\ ph a is ~ O.lpc 2 , 
whereas the effective cross sectional area for the self scattering case 
was just ~ 10~ 6 pc 2 . As such, if the scattering efficiency, ?7A,is 
large enough, this could be an important mechanism for growing 
black holes. 

From equations and H, 

we get for the mass accreted by 
hole B as black hole A moves from a radius a to b: 



AM E 



rj A 7YG 2 M A Vc(r) 4 
K <v 2 > 2 r' 



-dr 



(30) 



For a power law density profile (see equation |l6|) and using 
equation Il8, this gives us: 



AM B = A'^MaVaM-t) 





A' 



1 



0.856 ln(l + A 2 ) 



(31) 
(32) 
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3.2.1 The scattering efficiency 

Comparing equations 20 and 31, we can see straight away that, if 
the scattering efficiency, tja, were nearly 1 then nearly all of the 
central dark matter would be swept up as the holes spiral together. 
This is because the constant, A' is larger than A by a factor ~ c 2 . 
As such, the value of t\a is key to determining the magnitude of the 
accretion effect. 

We can obtain a crude estimate for t\a if we assume that all 
dark matter particles are scattered randomly. If this is the case then 
the probability, Pa, that a particle will be scattered into hole B's 
loss cone will be simply the fraction of phase space subtended by 
its loss cone. 

Thus we get, from equations 1,3, 17 and[l8|: 



Pa 



1 

4tt 



Q.A-KV 2 g(v)dv 



9GM 2 B (3-a) 



(33) 



(34) 



2TTC 2 p r™r<. 4 - a ) 

If Na is the number of times hole A crosses the same section 
of configuration space, then the fraction of particles scattered into 



B's loss cone will be: 



T) A = 1 - (1 - Pa) 



N a 



(35) 



Na is essentially the number of attempts hole A can make at 
scattering the same particles into hole B's loss cone. 

Thus, if the time taken for hole A to complete one orbit at a 
radius, r is given by T(r) — 2 Z r , , then Na is given by: 



dN A 



dt 



v c {rf 



Na = 



2nKM A r 2 p(r 

8A'(a)poro° 
M A (3-ay ( ° 



■dr 



(36) 



So that Na is now the number of times hole A crosses the 
same section of phase space as it moves from a radius a to a radius 
b. Thus, if a — b — Ar <C a, such that Pa(cl) — PaQ>) and if 
t]A < 1 then: 



Va(o) 



-N A P A 



1 — exp 



-36GA'(a)M 2 B (a- b) 
nc 2 M A (3-a)a 2 



(37) 



Notice that the scattering efficiency is maximised for Ma <C 
Mb - This is because the size of hole B's loss cone is proportional 
to M B (see equation 2) and so as Mb increases, the chance of hole 
A scattering matter into B's loss cone increases as M%. As Ma 
increases, the infall time for hole A decreases and so hole A will 
cross the same section of configuration space fewer times, reducing 
the scattering efficiency. 

The dependence of the scattering efficiency on a is seemingly 
very weak. However, A' is a strong function of a, particularly for 
small r (see figure |5| and equation 32). Since A' increases as a is 
reduced, we can expect that t\a will be similarly correlated. 

Notice also that r\A is independent of po and r . This is because 
increasing the central density normalisation por$, while keeping 
the black hole mass constant reduces the fraction of phase space 
subtended by the loss cone, and so Pa falls. However, the number 



10.0 
r(pc) 



Figure 4. The scattering efficiency t\a , as a function of radius for Ma 
M b = 1O 8 M ,« = 1. 



100.0 1000.0 




100.0 1000.0 



Figure 5. Mass accreted by hole B due to the scattering of dark matter 
particles by hole A. Graph (a) is for Ma = Mb = 1O 8 M0 . Graph (b) is for 
Ma = M B = 1O S M . Graph (c) is for M A = 10 2 M Q , M B = 10 8 M Q . 
Graph (d) is for M A = 10 2 M Q , M B = 10 6 M©. The solid lines are for 
a = 0.5, the dotted lines are for a = 1 and the dashed lines are for a = 1.5. 
The infall times, T, are marked in years for each case. 



of times hole A sweeps through the same section of phase space, 
Na, is subsequently increased. The two effects cancel out, so that 
the scattering efficiency depends only on the slope of the density 
profile and not its normalisation. 

Finally, notice that the scattering efficiency is a strong function 
of radius and will be maximised for small a. For most radii, the 
scattering efficiency is absolutely tiny. Figure |4] plots rjAifl) against 
a for Ma = Mb = 1O 8 M0, a = 1. As can be seen for almost 
all radii, as might be expected physically, the scattering efficiency 
is practically zero. Only when the holes are close does it become 
large enough to produce any significant amount of accretion. 



3.2.2 Results 

The numerical integration wasperformed using IDL, with the initial 
conditions set up as in table [l and with a numerical resolution of 
Ar = a — b — 0.1 parsecs. The results are displayed in figure |i[ 
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Figure 6. The change in the central dark matter density due to accretion 
onto two merging 10 8 Mq black holes. Graphs (a) and (c) show the central 
density before (solid line) and after (dotted line) the merger, for an initial 
density cusp of slope, a = 1 and a = 1.5 respectively. Graphs (b) and (d) 
show the central mass interior to r before (solid line) and after (dotted line) 
the merger, also for an initial density cusp of slope, a = 1 and a = 1.5 
respectively. 



As can be seen from equations 37, 31 and figure [j, accretion 
becomes significant only for small r. Reducing a, the slope of the 
cusp, pushes the radius at which accretion becomes significant to 
larger radii but does not significantly affect the total mass accreted. 

Figure 5(a) shows an extreme case -for Ma = Mb = 10 8 Mq. 
Here, some 10 5 Mq is accreted. While this amounts to only ~ 0.1% 
of the mass of hole B, it is sufficient to completely remove all of 
the dark matter within r c ~ 20 parsecs (see figure ^). Figure 5(b) 
shows a more likely scenario - for Ma = Mb = 1O 6 M0. For an 
NFW type inner slope of a=l or shallower, ~ IOOM0 is accreted. 
However for steeper slopes this rapidly falls to a mere O.1M , 
which is accreted only within the inner few parsecs. Figure 5(c) 
shows the case for Ma = 1OOM , M B = 1O 8 M . While only a 
few tens of solar masses are accreted, in the early stages of black 
hole growth, many low mass holes could have coalesced to form a 
supermassive central black hole. Miralda-Escude & Gould (2000) 
estimate that some 25,000 low mass black holes would have fallen 
into the central few parsecs of the Milky Way by the present day. 
If each of these scatters ~ 1OM of dark matter onto the central 
black hole then all of the dark matter within r c ~ 0.01 parsecs 
(~ 2.4 x 10 5 Mq) would be swept up. As with the case (a), this 
amounts to only ~ 0. 1% of the mass of hole B. Finally, Figure 5(d) 
shows the case for M A = 100M Q , M B = 1O 6 M . As with the 
case (c), many infalling black holes could lead to all of the dark 
matter within r c ~ 1 parsec or so being swept up. However, this 
amounts to only ~ 2OOOM , again only ~ 0.1% of the mass of 
hole B. 



However, even in the most extreme situations, the overall mass 
accreted is negligible, and usually less than that accreted in the 
steady state scenario. 

Mutual scattering is more effective, but only at very small 
radii. The main problem with this is that at radii of a few to a 
few tens of parsecs, the black holes are lik ely to form a stable 
hard binary (Milosavljevic & Merritt 2001). Dynamical friction 
then no longer applies since, as the binary hardens, much of the 
central matter is scattered away leaving little to disrupt the binary 
system. This mechanism, proposed by Quinlan (1996) and others, is 
effective on a scale of hundreds of parsecs. However, dense central 
clusters are found around massive black holes, so this process cannot 
be universal. Lauer et al. (2002) have recently identified several 
galaxies where central luminosity profiles are seen with cores of 
~ 20 parsecs. It seems likely, then, that well before any significant 
accretion could take place, much of the central matter would have 
already been removed by the formation of a hard supermassive black 
hole binary. This binary system would then decay, in the usual way, 
via gravitational radiation. 

The mutual scattering mechanism can, however, also be ap- 
plied to infalling massive satellites. A small bound system would 
survive right into the centre but be disrupted well before the forma- 
tion of a binary. It seems unlikely, however, that the satellite would 
remain bound as close as 10 parsecs to a supermassive black hole. 
The main point though is that, even if this were the case, the total 
mass accreted would still be a tiny fraction of the central black 
hole's starting mass and would affect a cold dark matter cusp only 
on the scale of a few parsecs. It seems that supermassive black holes 
cannot alter a central dark matter distribution on the kiloparsec scale 
required to align cold dark matter theory predictions with data from 
rotation curves. 

Finally, removing the constraint of spherical symmetry is un- 
likely to drastically alter our results. Although the model we present 
is rather crude, the magnitude of the accretion affect is found to be so 
small, that it is hard to imagine how a change in initial assumptions 
could lead to significant accretion. While it is true that particles 
lying on highly eccentric orbits are more likely to be accreted, the 
effect of strong dynamical friction as required in our model is likely 
to circularise the or bits of the stars and dark matter, reinforcing our 
initial assumptions (Polnarev & Rees 1994). 



This seems to contradict the value cited by ZHR of up to 
20-40% of the central black hole's mass being comprised of dark 
matter, however the reason for this discrepancy is clear. ZHR look 
at the effect of the diffusion of dark matter and stars into the loss 
cone via relaxation, whereas we look at refilling the loss cone by 
scattering from massive infalling bodies. As mentioned above, if 
the loss cone refilling rate is significant, then all of the dark matter 
within lkpc could be swept up leading to a significant fraction of a 
black hole's mass comprising of dark matter. In previous literature 
(see e.g. Lightman & Shapiro 1977) the diffusion rate of matter into 
the loss cone was calculated by making three main assumptions: 

1. The supermassive black holes are fed from shallow cores. 



4 DISCUSSION 

We have looked at some simple models which treat black holes in 
spherical systems as scattering bodies. A fraction of the scattered 
material is removed by accretion onto the scattering bodies, either 
through self scattering or mutual scattering as discussed above. 

Self scattering is most effective when the black holes are at 
their maximum separation, and for steep central density slopes. 



2. The stellar distribution is spherical or axisymmetric and angular 
momentum is conserved. 

3. Star-star relaxation is the dominant process for repopulating the 
loss cone. 

ZHR question these three assumptions in their paper and show 
that, the loss cone refilling rate could be significantly higher than 
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previously thought. They then go on to add that if there were no 
depletion of orbits for the dark matter, that is, if the loss cone were 
kept continually filled then some 20-40% of a supermassive black 
hole's mass could be comprised of dark matter. While this is true, 
it is somewhat misleading. 

The relaxation time for N particles in a collisionless system 
and assuming that eac h particle moves in the me an potential of all 
the others is given by (Binney & Tremaine 1987): 



N 



8 In TV \GM(r 



(38) 



Where M(r) is the mass interior to radius r, and TV is the 
number of particles. 

Now, given that the number of dark matter particles must be 
very very large, and since t re i ax oc TV, the dark matter-dark matter 
relaxation time will be very much longer than the age of the universe. 
However, a dark matter particle moving in the granular field of many 
stars will behave in the same way a particle of any mass moving 
through the same field. As such, the dark matter, even if cold and 
with very low self-interaction cross section, can be relaxed in the 
same way as the stars by the stars. 

The question is at what radius does the relaxation time for stars 
and dark matter fall below a Hubble time? If we simply model a 
galaxy as an exponential disk of stars, a central massive black hole 
and a power law dark matter profile, then the relaxation time is 
given by: 



.(»•) = 



TV, 



8 In TV, V G(Mbh + M b (r) + M dm (r) 



(39) 



Where Mbh is the black hole mass, Mb(r) is the mass distri- 
bution of the baryons and Md m is the mass distribution of the dark 
matter. We can approximate TV, ~ Mb(r) / Mq and so, for given 
mass profiles and a given central black hole mass, the relaxation 
time can be found as a function of radius. For a Milky Way type 
galaxy, the relaxation time is only less than a Hubble time within 2 
parsecs of the central black hole. The total matter out to this radius is 
not a significant fraction of the central black hole's start mass. ZHR 
argue that the relaxation time will be less than given in equation 
due to the relaxation of assumptions 1-3 above. They present some 
plausible arguments that could reduce the relaxation time and then 
state that, should it be sufficiently low that the loss cone can be kept 
refilled out to the black hole's sphere of influence, then 20-40% of 
a central supermassive black hole's mass could comprise of dark 
matter. It is not clear, however, whether the relaxation time really 
could be low enough to produce significant loss cone refilling. Even 
if it were, the radius to which the effect can be considered important 
must be less than the sphere of influence of the black hole, and so 
cannot be important for producing dark matter cores on the scale 
of a kiloparsec or so - the kind of core which would be required to 
reconcile measurements from rotation curves with cold dark matter 
theory. Finally, they do not consider the effect of mergers on the 
central matter profile. As stated above, a merging binary black hole 
can eject much of the central matter from a galaxy on the scale of 
10-100 parsecs or so. Even if the relaxation time at 100 parsecs 
were quite low, it is likely that the matter in this region would be 
removed by merging black holes before it would have the chance 
to be accreted - either by diffusion due to the coarse grained stel- 
lar potential, as suggested by ZHR, or by scattering from merging 
supermassive black holes as suggested here. 



5 CONCLUSIONS 

We have shown that supermassive black holes cannot alter a power 
law density cusp via accretion, whether during mergers or in the 
steady state. This represents another vital piece in the debate sur- 
rounding the nature of dark matter. If winds and bars also prove to 
be ineffective at altering a primordial dark matter cusp, we will be 
drawn inexorably towards the conclusion that dark matter is more 
complex than the simple cold particle envisaged so far. 
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7 APPENDIX I - THE VELOCITY DISPERSION IN A 
POWER LAW DENSITY PROFILE 

For sphericall y symmetric, isotropic, systems we have (B inney & 
Tremaine 1987): 



1 d(p < v- >) 
P 



dr 



Vc 

r 



(40) 



If the density, p(r), is given by equation |l^ then v\ is given by 
equation 17 and equation ^0j becomes: 



_i 2 d<ir > 

-OCT < V > +- 



dr 



-4nGp r^ rl _ a 



(41) 



< v > 



2-KGpprp 
—1-kG p^r^x In x 



1 < a < 3 



(42) 



Since the power law density is describing the inner parts of 
the halo, we require r < ro always and so x < 1 . Thus In x will 
always be negative such that < v 2 > is still positive for a = 1. 
Thus substituting everywhere for v\ (see equation 17) we get: 



< v > 



< v > = 



2(o-l) 

v c 1 < a < 3 

—vl lax 

2 , 
v c a — 1 



(43) 
(44) 
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